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The universal behaviour of two-dimensional loop models can change dramatically when loops 
are allowed to cross. We study models with crossings both analytically and with extensive 
Monte Carlo simulations. Our main focus (the 'completely packed loop model with crossings') 
is a simple generalisation of well-known models which shows an interesting phase diagram with 
continuous phase transitions of a new kind. These separate the unusual 'Goldstone' phase observed 
previously from phases with short loops. Using mappings to Z2 lattice gauge theory, we show that 
the continuum description of the model is a replica limit of the sigma model on real projective 
space (RP n_1 ). This field theory sustains Z2 point defects which proliferate at the transition. 
In addition to studying the new critical points, we characterise the universal properties of the 
Goldstone phase in detail, comparing renormalisation group (RG) calculations with numerical data 
on systems of linear size up to L — 10 6 at loop fugacity n = 1. (Very large sizes are necessary 
because of the logarithmic form of correlation functions and other observables.) The model is 
relevant to polymers on the verge of collapse, and a particular point in parameter space maps to 
self-avoiding trails at their O-point; we use the RG treatment of a perturbed sigma model to resolve 
some perplexing features in the previous literature on trails. Finally, one of the phase transitions 
considered here is a close analogue of those in disordered electronic systems — specifically, Anderson 
metal-insulator transitions — and provides a simpler context in which to study the properties of 
these poorly-understood (central-charge-zero) critical points. 



I. INTRODUCTION 

It is tempting to think that two-dimensional criti- 
cal phenomena are completely classified and understood, 
thanks to conformal field theory and other exact tech- 
niques, but this is far from true. One class of problems 
which is still mysterious contains classical loop models 
and models for polymers, together with models for non- 
interacting fermions subject to disorder. These systems 
are tied together by field theory descriptions with con- 
tinuous replica-like symmetries (or alternatively global 
supersymmetries). Examples include de Gennes' map- 
ping of polymers to the O(N) model in the limit N — >• 0; 
the various replica sigma models for the integer quantum 
Hall transition and other localisation problems; and those 
describing cluster boundaries in percolation and similar 
soups of loops [ThT5] . 

Classical loop models yield the simplest examples of 
this class of problems, but even they are not yet fully un- 
derstood. The best-studied examples are those in which 
the loops are forbidden from crossing: for these a great 
deal is known from conformal field theory, height model 
mappings, exact solutions, Schramm-Loewner Evolution, 
and numerical simulations [141 115] . But when we move 
away from these models the analytical techniques often 
cease to apply, and we may encounter new types of crit- 
ical phenomena requiring new theoretical tools. 

This paper considers some of the simplest two- 
dimensional loop models with crossings. These reveal 
new universality classes of, and new mechanisms for, clas- 
sical critical behaviour. They also provide natural models 
for polymers and for deterministic motion in a random 



environment [16-18] which have been intensely studied 
but whose phase diagrams and continuum descriptions 
have in general not been understood. Finally, they shed 
light on phenomena that are important more generally 
for criticality in replica or supersymmetric sigma mod- 
els - in particular, the role of gauge symmetries and 
topological point defects. The latter have recently been 
shown also to be important for two-dimensional Ander- 
son metal- insulator transitions [121 [2U] . We will return 
to the analogy between loop models and localisation at 
the end of this introduction. 

A key result of previous work on loops with crossings 
is the existence of an unusual critical phase which is ab- 
sent for non-crossing loops [TUl [TT ] [TCI \T7\ I2TH25] . It was 
argued by Jacobsen, Read and Saleur [TUl E] that this 
corresponds to the Goldstone phase of the 0(n) sigma 
model, where n is the fugacity for loops. The phase ex- 
ists for n < 2; to make sense of this regime requires 
a replica-like limit or a supersymmetric formulation of 
the field theory. Characteristic features of the Goldstone 
phase had previously been found in computational stud- 
ies of polymers and deterministic walks in a random en- 
vironment [16j [17] , as well as in an integrable loop model 
[2TM23] . The phase appears quite generically when non- 
crossing loop models in the so-called 'dense' regime are 
perturbed by the addition of crossings, which in an ap- 
propriate field theory corresponds to a breaking of sym- 
metry [TUHH]- 

Here we examine a more general class of loop models 
with crossings. These show new continuous phase tran- 
sitions separating the Goldstone phase from non-critical 
phases with short loops. We construct field theories for 
the models, and pin down the universal behaviour (both 
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FIG. 1. A configuration of the compietely packed loop model 
with crossings (CPLC) on a 10 x 10 lattice with periodic 
boundary conditions. Each loop has been given a different 
colour - these colours are merely a guide to the eye and not 
part of the configuration. 



in the Goldstone phase and at the new critical points) 
using analytic calculations and extensive Monte Carlo 
simulations. Finally, we give a field theoretic treatment 
of the closely related problem of the interacting self- 
avoiding trail model for a polymer. 

The models we study are a 'completely-packed loop 
model with crossings' (CPLC) on the square lattice - 
Fig. [T] shows a configuration — and an 'incompletely 
packed loop model with crossings' (IPLC) in which 
loops are related to cluster boundaries. The param- 
eter space for the CPLC contains various previously- 
studied models as special cases: in particular, the stan- 
dard completely-packed loop model without crossings and 
models with crossings encountered in various contexts 
(including those mentioned above). The CPLC and IPLC 
are expected to show the same universal behaviour — our 
numerics are restricted to the CPLC, but the IPLC pro- 
vides a simpler context in which to describe the main 
theoretical ideas. 

It is easy to argue that the phase transitions in the 
CPLC and IPLC cannot be described by the 0(n) model. 
Instead, the general description of the models requires 
couping the O(n) spin S to a Z2 gauge field, or equiva- 
lently identifying S with —S to obtain a nematic order 
parameter. This leads to a sigma model on real projec- 
tive space, RP™ _1 , in which Z2 point vortices play an 
important role. Vortices are suppressed in the Goldstone 
phase — meaning that the 0(n) model is a viable de- 
scription there — but proliferate at the phase transition 



into the short loop phase. (This picture is appropriate 
for the regime < n < 2.) 

In general, the introduction of crossings leads the stan- 
dard exact techniques used for non-crossing loop models 
to fail, so for the critical points we are restricted to nu- 
merics and approximate RG treatments j20l . However 
the Goldstone phase can be understood fully, since it is 
characterised by marginal flow to a weak-coupling fixed 
point [TUJ [Tl] . This leads to logarithms (e.g. correlation 
functions decaying with a universal power of the loga- 
rithm of distance), so in order to confirm the analytical 
predictions, very large sizes are required for the numerics 
(comparable to the largest system sizes simulated in any 
statistical mechanics problem). These are possible at fu- 
gacity n = 1 thanks to special features of the problem 
there, and our simulations are restricted to this value. 

Another feature of the CPLC at n = 1 is that while 
each configuration is a soup of many loops, the model 
permits a mapping to a model for a single loop with local 
interactions. At a certain point in parameter space, this 
is the well-studied 'interacting self-avoiding trail' (ISAT) 
model for a polymer at its collapse, or 6, point [TBI |2"1] . 
Collapse transitions for polymers in two dimensions are 
a mysterious subject which deserves clarification (see for 
instance the case of the missing Flory exponents [25j). 
Here we show that the ISAT can be understood com- 
pletely from field theory, explaining for example the in- 
teresting phase diagram found numerically in Ref. [18) . 
To do this we perturb the sigma model which describes 
the Goldstone phase. Surprisingly, the point of the 
ISAT turns out to be an infinite order multicritical point: 
despite the simplicity and naturalness of this model, it is 
highly fine-tuned from the point of view of general poly- 
mer models. This implies that the critical exponents for 
the generic point polymer (with crossings) are still 
unknown. 

There is a close relationship between loop models at 
loop fugacity n — 1 and disordered fermion problems |26l - 
I32j . Supersymmetry and replica- like limits, crucial in the 
latter for averaging over disorder, appear in the former 
as tools allowing geometrical correlation functions (such 
as the probability that two distant sites lie on the same 
loop) to be expressed. Both types of problem exhibit 
critical points of central charge zero, described by loga- 
rithmic CFTs [33]. The loop models are a good place to 
study such critical points since they are more tractable, 
both analytically and computationally, than disordered 
fermion problems. 

For completely-packed loops without crossings, there is 
in fact an exact mapping |26l 12 7j to a network model for 
Anderson localization in symmetry class C [3H 155] . 
However the analogy is more general. Recent work by 
Fu and Kane [5D] demonstrates that the metal-insulator 
transition in the symplectic symmetry class is driven by 
proliferation of Z2 vortices: this transition is thus in re- 
markably close analogy with those in the loop models dis- 
cussed here, though the appropriate sigma model is dif- 
ferent. In the localisation language, the Goldstone phase 
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corresponds to a metallic phase, and the two short loop 
phases — which are distinguished from each other by the 
presence or absence of a loop encircling the boundary — 
to topological and trivial insulating phases. 

The CPLC at n = 1 can in fact be obtained as a 'classi- 
cal' limit of a network model in which a Kramers doublet 
propagates on every edge: the above similarities show 
that this classical limit captures a surprising number of 
the qualitative features of the phase diagram for the sym- 
plectic class. 

In both the loop model and the localisation problem 
the vortex fugacity (including its sign) plays an impor- 
tant role. Ref. [20] introduced an approximate RG treat- 
ment of this fugacity, and in Sec. V C| we apply this to the 
loop models. We note that vortices (this time Z vortices) 
have also been shown to be responsible for Anderson lo- 
calisation in the chiral symmetry classes, and a detailed 
treatment has been given by Konig et al. in Ref. 19J. 

The organisation of the paper is as follows. In the 
next section we introduce the models we will study and 
their phase diagrams. In Sec. Ill we map them to lattice 
gauge theories and RP™" 1 sigma models, paying atten- 
tion to the role of topological defects and the relation 
to the field theory for loops without crossings. We also 
discuss the different (CP™ -1 ) sigma model which applies 
on the boundaries of the phase diagram for the CPLC. 
Sections [TV] and [V] apply Monte Carlo and RG calcula- 
tions in the field theory to the Goldstone phase and the 
critical points respectively. Our numerical methods are 
described in more detail in Sec. EE] Sec. IVHI tackles the 
polymer collapse problem. Finally, Sec. |VIII| discusses 
directions for future work. 



II. DEFINITIONS OF MODELS 
A. Completely-packed loops with crossings 

Configurations of the completely packed loop model 
with crossings (CPLC) are generated by resolving each 
node of the square lattice in one of the three possible ways 
shown in Fig. [2} Fig. [TJ shows an example on a small lat- 
tice. Note that overcrossings are not distinguished from 
undercrossings — the configuration at a node is defined 
solely by the way its four links are paired up. 

Each of the three possible pairings at a node is assigned 
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FIG. 2. The three configurations of a node and associated 
Boltzmann weights. (In the leftmost configuration, the upper 
and lower links lie on the same loop.) The weights q and 1 — q 
are exchanged on the two sublattices of the square lattice. 



a weight, as shown in Fig. [2] with the weight of a crossing 
being p. The factors q and 1 — q are staggered (swapped) 
on the two sublattices of the square lattice, so that the 
states of the system for extreme values of the parameters 
are as shown in Fig. [3] The Bolt zmann weight for a 
configuration is given by the product of the node factors, 
together with a fugacity n for the loops. Let N p , N q and 
Ni- q denote the numbers of nodes where the pairing with 
weight p, (1 — p)q or (1 — p)(l — q) is chosen. Then the 
product of node weights in a configuration C is 



W C =P N * [(l-p)q] N <[(l-p )(!-'/> 
and the partition function is 

Z = J2 lo ° pS W / c- 
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The parameter space of this model includes various 
previously-investigated models. On the line p — we 
have the completely-packed loop model without cross- 
ings, which has been intensely studied and which may 
be mapped to the n 2 -state Potts model via the Fortuin- 
Kasteleyn representation of the latter. The model on 
the line q = 1/2 was related to the Goldstone phase of 
the 0(n) sigma model in Ref. [TO], and points on this 
line have been studied in various contexts. For a given 
value of n the point q = 1/2, p = (2 — n)/(10 — n) is 
known as the Brauer loop model [2"TH2"3l and is integrable; 
this model was related to a supersymmetric spin chain in 
Ref. [21]. When the parameters in the CPLC are such 
that all configurations given equal weight — i.e. when 
n = 1, q = 1/2 and p — 1/3 — it is equivalent to a 
standard model for polymers at their O (collapse) point 



[TBI fT8l [24], which we will discuss further in Sec. VII 
On the lines g = or g = 1 — the left and right bound 
aries of the phase diagram Fig. [3]— the CPLC reduces to 



the 'Manhattan' lattice loop model discussed in [271136] . 
Loop models with crossings at n = 1 have also appeared 
in the study of Lorentz lattice gases, i.e. determinis- 
tic motion in a random environment [171 137j . Finally, 
Ref. [38] discusses a model similar to the CPLC in which 
q and 1 — q are not staggered, and uses it (at n = 2) to 
analyse the phase diagrams of vertex models. 

A trivial but important fact about the CPLC is that 
the nodes become completely independent of each other 
when n = 1. The weights p, (1 — p)q and (1 — p)(l — q) 
are then the probabilities of the various node configu- 
rations, and the partition function Z is equal to unity 
— from which it follows, by the finite size scaling of the 
free energy, that any critical points must have central 
charge c = 0. The model with n = 1 is thus analogous to 
percolation, which can also be formulated in terms of un- 
corrected random variables. In the absence of crossings 
the n — 1 model is in fact equivalent to bond percolation 
on a dual lattice, with loops surrounding cluster bound- 
aries; however when crossings are allowed the universal 
behaviour is no longer that of percolation. 

Our simulations will be restricted to the case n = 1, 
which is the most interesting and the best suited to 
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FIG. 3. Phase diagram obtained numerically for the CPLC 
at n = 1. The horizontal axis is labelled by q, defined by (q — 
1/2) = (q — 1/2) (1 — p). The larger (red) dots on the critical 
line indicate the values of p at which we have analysed the 
critical behaviour in detail. Also shown are the configurations 
obtaining on a small finite lattice at p = 1, at p = 0, q = 0, 
and at p = 0, q = 1. (The point p — 0, q = 1/2 is the 
percolation critical point.) 



Monte Carlo, but most analytic results will apply to 
< n < 2. (We will discuss the case n = 2, which 
shows more conventional critical behaviour, elsewhere.) 
The phase diagram obtained numerically at n = 1 is 
shown in Fig. [3| We expect it to be qualitatively similar 
for < n < 2, with the Goldstone phase swallowing up 
more and more of the parameter space as n -> 0. We 
now summarize its main features. 

Short loop phases. The configurations at p = 0, 
q = and at p = 0, q — 1 provide caricatures of the 
two 'short loop' phases. For given boundary conditions, 
these are distinguished from each other by the presence 
or absence of a long loop running along the boundary, 
as shown in Fig. |] In a sigma model description the 
short loop phases are massive (disordered) phases. In 
the analogy with Anderson localization mentioned in the 
introduction, they correspond to insulating phases, and 
the boundary loops correspond to the edge states present 
in a topological insulator. 

Goldstone phase. In the Goldstone phase, so-called 
because the continuum description is a sigma model 
which flows to weak coupling in the infra-red [TU], the 



loops are 'almost' Brownian. However, interactions be- 
tween Goldstone modes in the sigma model are only 
marginally irrelevant, leading to universal logarithmic 
forms for correlators and other observables which we cal- 



culate in Sec. IV In the Anderson localization analogy, 
this would be a metallic phase. 

Critical lines. The lines separating the Goldstone 
phase from the short-loop phases show a new universal- 
ity class of critical behaviour. This is associated with 
the order-disorder transition of the RP™ -1 sigma model, 
which exists only in the replica limit n < 2 and is 
driven by proliferatio n of Z? vortex defects associated 



with TTl 
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V C I . Numerically, the critical 
loops have df = 1.909(1) at n = 1, i.e. they are slightly 
less compact than Brownian paths, and the transition 
has a large correlation length exponent v = 2.745(19) 
(Sec.fv). 

Critical point at p = 0. The critical point of the 
loop model without crossings (at p = 0, q = 1/2) is well- 
studied and corresponds to the so-called dense phase of 
the 0(n) loop modefj] , or to SLE K for k > 4. At n = 1 
the model maps to critical percolation: the loops have 
the statistics of percolation cluster boundaries, with a 
fractal dimension dj GIC = 7/4, and the correlation length 
exponent of the transition is i/ pcrc = 4/3. These critical 
exponents also yield exponents in an Anderson transition 
(the spin quantum Hall transition) via an exact mapping 
[26H2H1I39]. 

Phase diagram boundaries. Everywhere on the 
boundary of the phase diagram Fig. [3] — i.e. whenever 
one of the node weights vanishes — loops can be consis- 
tently oriented by assigning a fixed (configuration inde- 
pendent) orientation to each link of the lattice. The nec- 
essary choice of orientations differs on each of the three 
pieces of boundary. They are those of the 'L' lattice 
on the line p — 0, and of the Manhattan lattice on the 
lines q = and q = 1. (These lattices are depicted in 
Ref. [27].) The fact that the loops automatically come 
with an orientation means that the continuum descrip- 
tions have a higher symmetry [21 HO] , as will be discussed 
in Sec. |III F| and the appropriate field theory is a sigma 
model on CP"" 1 rather than on RP"" 1 . The CP" -1 de- 
scription implies that the lines q — and q = 1 (the 
Manhattan lattice loop model) are always in the short 
loop phase, but with a typical loop length that diverges 
exponentially as p — >• 1 (Sec. |IIIF ). This is in agreement 
with previous expectations 27, 13(3] . but is not obvious 
from the numerical phase diagram since the critical lines 
closely approach the lines q — 0, 1 for p close to one. 

Relation to polymers. Configurations in the CPLC 
are soups of many loops. However when n = 1 the CPLC 
has a simple relation with the self-avoiding trail model 
for a single polymer [161 118] . The polymer corresponds 



1 Note that the standard terminology 'O(n) loop model' is poten- 
tially misleading, as when crossings are forbidden the 0(n) model 
is not the appropriate field theory for the dense phase [9] I10| . 
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FIG. 4. Part of a configuration in the IPLC. In this model 
the loops (thick/red lines) are cluster boundaries. Universal 
behaviour in the IPLC is expected to coincide with that in 
the CPLC. 



to a single marked loop in the soup of loops; so long as 
n = 1, 'integrating out' the configurations of the other 
loops leads to a local Boltzmann weight for the marked 
one. Adding the interactions which are natural in the 
polymer language takes us out of the parameter space of 
Fig. |J but the sigma model description can be extended 
to cover this case by including appropriate symmetry- 
breaking terms (Sec. VII). 



probability, the partition function for the IPLC is 

z = 2^ a q i 1 - q) > 

configs 



(3) 



with a — 1/3. The parameter q here will play a similar 
role to the parameter q in the CPLC. 

The above partition function corresponds to a loop fu- 
gacity n = 1. We will wish to generalize it to arbitrary 
n. We may also vary a, and introduce a fugacity x for 
the total length of loops: 



Z 



c 



(1-q) 



W„length no. loops 



(4) 



B. Incompletely-packed loops with crossings 



The precise values of a and x will not be important in 
what follows. When n = 1, a length fugacity x distinct 
from one corresponds to an Ising interaction of strength 
J between colours of adjacent squares, with x = e~ 2J , 
and varying a introduces a four-site interaction. These 
interactions have no effect on the universal behaviour so 
long as they are weak. 

Note that if we take our lattice to have the topology 
of the disk, regarding the region outside the boundary 
as white for the purposes of drawing cluster boundaries, 
the IPLC shares with the CPLC the feature of having an 
edge loop at q = 1 but not at q = 0. There are clearly 
stable short-loop phases at small q and at q close to one, 
and field theory arguments lead us to expect a stable 
Goldstone phase near q = 1/2. The conjectured phase 
diagram, shown in Fig. [5] is similar to a slice through the 
phase diagram of the CPLC at some nonzero value of p. 
(If we had forbidden crossings, we would have obtained 
a phase diagram similar to the line p = in the CPLC.) 



For pedagogical reasons, it will be useful to introduce 
and discuss a second loop model before returning to the 
CPLC. Loops in the new model will no longer be com- 
pletely packed, but nevertheless the universal proper- 
ties will be the same. We refer to this model as the 
incompletely-packed loop model with crossings, or IPLC. 

To generate a configuration in the IPLC, we first colour 
the plaquettes of the square lattice grey or white, giv- 
ing a site percolation configuration on the square lattice 
formed by the plaquettes. The loops in the IPLC are then 
cluster boundaries, as shown in Fig. |4j However the loop 
configuration is not uniquely determined by the cluster 
configuration: for each 'doubly visited' node, where two 
cluster boundaries meet, we must choose how to connect 
them up. Allowing crossings, the three possible pairings 
are again those of Fig. |] (but unlike in the CPLC we do 
not assign different weights to the different pairings). 

The simplest choice for the Boltzmann weight is to 
give each percolation configuration the standard perco- 
lation probability q B (l — q) w , where B and W are the 
numbers of black and white faces. A given percolation 
corresponds to 3 N loop configurations, where N is the 
number of doubly-visited nodes. Assigning them equal 



III. LATTICE FIELD THEORIES 
A. Lattice field theory for IPLC 

We begin with the IPLC, which permits a simple map- 
ping to a Z 2 lattice gauge theory coupled to matter fields. 
We will first write down this theory and then show that a 
graphical expansion similar to the high temperature ex- 
pansion of the Ising model or the Nienhuis 0{n) model 
[HI p4"2"] provides the connection to the loop model. 

The required lattice field theory includes matter fields, 
which are real n-component vectors living on the sites i 
of the square lattice, 



S 



(5) 



Short loops 1 



Goldstone 



Short loops 2 



FIG. 5. Schematic phase diagram for the incompletely-packed 
loop model with crossings as a function of q (for fixed n ~ 1, 
a ~ 1/3, x ~ 1). 
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and gauge fields cry 
tion function is: 



Z = TSrJj((l-q)+q J] **) U 



±1 living on the links. The parti- 



' X&ij Si-Sj 



'Tr' denotes sums and integrals over all degrees of free- 
dom (cr and S), normalized so Tr 1 = 1, and F denotes a 
face (square) of the lattice. The Z2 gauge symmetry of 
this model is 



(for x% 



±1) 



(7) 



Si -> XiSn a ij XiXjCij 

Above we have written the Boltzmann weight for the 
gauge field in a form suitable for the graphical expansion. 
Later we will rewrite it in a more conventional form. 

The graphical expansion begins by expanding out the 
product over faces F and representing the terms by a 
simple graphical rule (a face is coloured grey if the q 
term is chosen, and white if the 1 — q term is chosen). 
This generates percolation configurations V, 




x&ijSi. S- 



(8) 

where &P denotes the set of links lying on percolation 
cluster boundaries. 

Next we expand out the product over links. Only one 
term in this expansion survives after summing over cr^-, 

namely that in which the factors of Si-Sj lie on the cluster 
boundaries: 

Tr Yl S i- S j- ( 9 ) 

V (ij)edV 

'Length' refers to the total length of cluster boundaries. 
We have written the spin index a explicitly in the inner 
product Si.Sj to emphasize that each link on a cluster 
boundary now carries an index a — 1, . . . n. Finally we 
perform the remaining integrals over the vectors S, using 

Tr = S a \ 



TV Q a Q c Q a 

-LI LJA LJ A iJA iJ A 



(S ab S cd + S ac S bd + 5 ad S bc ) . (10) 



The three terms in the second formula correspond to the 
three ways of connecting up the links in pairs at a node 
for which all four links are in cfP . Expanding out all 
such brackets gives 3^ terms, each associated with a loop 
configuration C. For a given C, each loop comes with a 
product of delta functions forcing the indices a to be 
equal for all links on that loop. We may think of each 
loop as carrying a 'colour' index ranging from l,...,n. 
Summing over the possible colours for each loop yields a 
fugacity n: 



c 



a' 



V(l-q) 1 



(11) 



The parameter a is n/(n + 2) as a consequence of Eq. 10 
It can be varied by exchanging the hard constraint S 2 = 
n for a potential for S 2 . 



B. Correlation functions and the replica trick 

We have seen that a graphical expansion of the lat- 
tice gauge theory ^ generates the configurations of the 
loop model with the right weights. The correspondence 
also extends trivially to correlation functions. Switching 
temporarily to a continuum notation, the basic ones are 
the watermelon correlation functions Gf.(x, y), which give 
the probability that k distinct strands of loop connect x 
and y. For example, on the lattice Gi can be taken to be 
the probability that two links lie on the same loop, and 
G4 the probability that two nodes are connected by four 
strands. 

In the continuum, the correlator is the two-point 
function of the 'ft- leg' operator Ok- 



O k (x) oc S\x)S 2 {x)...S k (x). 



(12) 



Gfc vanishes for odd k, either by gauge invariance or 
equivalently because two regions cannot be joined by an 
odd number of strands (for example if two sites lie on 
the same loop they are joined by two strands). We can 
visualize the above operator, when inserted into a corre- 
lation function, as emitting k strands with colour indices 
ranging from 1 to k (Fig [6]). 

We see that in order to write down the correlator G k 
we require the number n of spin components to be at 
least k. This presents a problem if the model we wish to 
study has a loop fugacity n < k: for the model of most 
interest to us with n = 1, it does not allow us to write 
down any of the above correlation functions! Fortunately, 
this problem can be resolved in two standard ways. 

Firstly, we can treat n — > 1 as a replica-like limit, so 
that 



G k = lim (O k (x)O k (y)) 

n— J-l 



(13) 



This trick may also be used for other values of the loop 
fugacity, integer or noninteger. For integer n there is 
a more rigorous alternative, which is to use supersym- 
metry [3J [5] . This allows us to increase the number of 
components of S without increasing the loop fugacity, 




FIG. 6. The watermelon correlator Gi(x, y). 
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by making some of the components fermionic. For our 
purposes the two approaches are equivalent, so for pre- 
sentational simplicity we will use the replica language. 
The required supersymmetric construction is explained 
in Refs. [TQl H2]: in outline, S is replaced by a vector $ 
with both bosonic and fermionic components, 



S 2 + 2ff.£ = n, 



(14) 



where S has n+2m components and ff, £ each have m an- 
ticommuting components. The loop fugacity determines 
n, but m is arbitrary - supersymmetry ensures that the 
partition function and its loop representation do not de- 
pend on it. Thus m may be chosen large enough for any 
desired correlator to be constructed. 



C. Z2 vortices and Z2 fluxes 

The local gauge-invariant degree of freedom in the lat- 
tice gauge theory ^ is a ncmatic vector, obtained by 
identifying S with —S. It can be encoded in a real sym- 
metric matrix, 



= S a S b - S ab 



tr Q = 0, 



(15) 



which will be the relevant degree of freedom on long 
length scales. Q lives on real projective space, RP" -1 : 
since this manifold has nontrivial fundamental group, the 
Q configuration can have vortex defects which we now 
discuss. 

RP 1 is equivalent to the circle, so at the special value 
n = 2 vortices are standard XY vortices and are char- 
acterized by an integer topological charge. However 



\ 



FIG. 7. Plaquettes with gauge flux H ff = ^1 (shaded in 
pink) are endpoints of strings of links with a — — 1 (marked 
in bold/blue), across which 5* changes sign. In terms of the 
nematic order parameter, which is obtained by identifying S 
with — S and which lives on RP n_1 , these plaquettes are vor- 
tices. (This is of course a caricature, neglecting fluctuations.) 



7Ti(RP™~ ) = Z2 when n > 2, so the vortex charge is in 
general defined only modulo two. In the replica limit — 
which requires analytic continuation of formulae defined 
for arbitrarily large n to n < 2 — the vortices should 
again be viewed as Z2 vortices. (This may be clearer in 
the supersymmetric formulation of the n = 1 loop model, 
Eq. |14| where the bosonic part of the superspin lives on 
IRP m for m > 1, which has fundamental group Z2.) 

For a more concrete picture we return to S and a. Let 
a/ = +1 on all links I, except for a semi-infinite string of 
parallel links ending at a plaquette F (see Fig. [7| . The 
flux Y[ v is then — 1 only on F. The spin configurations 
S which maximize the Boltzmann weight vary smoothly 
except at the string, across which S changes sign. This 
indicates the presence of a Z2 vortex located at F. (This 
connection between vortices in nematics and Z2 fluxes is 
standard [13].) 

Since vortices are associated with plaquettes of non- 
trivial flux, we can assign a vortex number to each pla- 
quette which is +1 if Y\a — —1 and if Y\a = 1. The 
gauge field part of the Boltzmann weight can then be 
written in terms of the number N v of vortices: 



J] (i-q) + q II a ^ =(i-2<0 



N, 



(16) 



(ij)eF 



We see that the factor (1 — 2q) is simply a fugacity for 
vortices, and that the exchange q •<-> 1 — q corresponds 
to changing the sign of the vortex fugacity. 

This sign distinguishes the two short-loop phases in 
Fig. [5] from each other. In Ref. [20], Z2 vortices play an 
analogous role in a sigma model for localisation, with the 
sign of the vortex fugacity distinguishing two insulating 
phases. (The vortex fugacity is also important for the 
critical behaviour [20] — see Sec. VC| ) 

Let us rewrite the Boltzmann weight for a in the con- 
ventional form for Z2 gauge theory. In the absence of a 
boundary, 



Z oc Tr exp k 



n 



Y[ (l + xaijSi.Sj 



f (ij)eF ' (ij) 
where the gauge field stiffness is 

1, _J_ 

k = - m r. 

2 |l-2q| 



(17) 



(18) 



In the presence of a boundary, denoted d, the Boltzmann 
weight acquires an additional term when q > 1/2: 



n < 

(ij)ed 



(19) 



This term effects the sign change in the vortex fugacity. 
It is equal to (-) a,s '™e ! j where iV s t r ings i s the number of 
a = — 1 strings which terminate on the boundary. Since 
this number is equal to the number of vortices in the 
interior modulo 2, (-) JVs,ri °s s = (~) Nv . 



We see that the sign of the vortex fugacity does not 
effect bulk properties. Instead it corresponds to the pres- 
ence or absence of an edge loop. 

Finally, consider the point q = 1/2. The vortex fugac- 
ity vanishes here (Eq. 16); however, the universal prop- 
erties of this point do not differ from those in the rest 
of the Goldstone phase (Fig. [5]) . This is because vortices 
are anyway RG irrelevant in that phase. We discuss this 
in the next section in terms of a sigma model for Q. 

The suppression of vortices (either microscopically or 
in the infrared) means that in the Goldstone phase this 
sigma model has a correspondence with the simpler 0(n) 
sigma model. In the IPLC at q = 1/2 this correspondence 
holds microscopically, since the gauge field stiffness di- 
verges there (18). This enforces \\a = 1 for every face 
F, giving (Tij — XiXj ( so l° n g as the lattice lives on a 
simply-connected manifolcQ 



Zcx^Tr^n 

fx} <«) 



l+x(xiSi).(XjSj) 



(20) 



Changing variables to S' — x$ eliminates the gauge de- 
grees of freedom from the Boltzmann weight, leaving 
a lattice 0(n) model. (Non-gauge-invariant correlators 
pick up factors of x> ensuring that they vanish^] on sum- 
ming over x-) I n t ne regime we consider, i.e. at suffi- 
ciently large x, this 0{n) model is expected to be de- 
scribed by the 0(n) sigma model in its Goldstone phase. 



D. Continuum description 

The naive continuum description of the lattice field 
theory ([6| is a sigma model for Q, 



C = ^tr(VQ) 2 



flows to large values under coarse-graining, and vortices 
are an irrelevant perturbation. 

Non-singular (i.e. vortex-free) configurations of Q are 
equivalent to non-singular configurations of S (on a man- 
ifold of trivial topology, and up to a global sign ambigu- 
ity). Thus for a perturbative treatment at large K, the 
RP™ _1 sigma model can be replaced with the more famil- 
iar 0(n) sigma model, argued previously [TO] to apply to 
the CPLC at q = 1/2: 



£=f(V5) 2 , 



S 2 = 1. 



(22) 



The perturbative beta function for K changes sign at 
n = 2 {MM- 



dK 
dlni 



2tt 



1 



2irK 



(23) 



(21) 



When n > 2, the stiffness flows to zero under RG and 
the sigma model has only a disordered phase; thus the 
loop model is not expected to be critical. At n = 2 the 
sigma model is the XY model, and we have in addition 
the quasi-long-range-ordered phase in which K does not 
flow. By the Mermin Wagner theorem, these are the only 
possibilities when n > 2. 

However in the replica limit the Mermin Wagner theo- 
rem does not apply [10] , and Eq. 23 shows that for n < 2 
the stiffness K flows to infinity in the infrared. This is the 
Goldstone phase. At the infra-red fixed point the n — 1 
Goldstone modes are free fields, so the central charge is 
c = n-l [TU]. 

We now discuss the extraction of a continuum descrip- 
tion for the CPLC, the appearance of a nontrivial vortex 
fugacity in that model, and the extra symmetry on the 
boundaries of the phase diagram. 



together with the constraint on Q following from its defi- 
nition in terms of S. To conform with convention we use 
the normalisation S 2 = 1 in the continuum (rather than 
Eq Is}; then Q ab = S a S b - ±5 ab . 

The fugacity for vortices is hidden in the ultraviolet 
regularisation of the sigma model (21). We will restore 



this parameter explicitly when we consider RG in the 
vicinity of the critical point (Sec. VC), where vortices 
are crucial. 

However, since the classical free energy of a pair of 
vortex defects is proportional to the stiffness K, they are 
suppressed at large K (as in the XY model at large stiff- 
ness). In the Goldstone phase, which we now discuss, K 



E. Lattice field theory for the CPLC 

The CPLC can again be mapped to a lattice spin model 
with a l-i gauge symmetry, though with a less conven- 
tional form, and as long as we are in the interior of the 
phase diagram (Fig. [3]) the continuum description is again 
the RP n_1 model. The following construction is similar 
to that described in Ref. [T3] so we give only an outline. 

We again introduce fixed-length spins Si (Sf = n) but 
for the CPLC they live on the links I of the square lattice. 
The Boltzmann weight involves four-spin interactions be- 
tween the spins surrounding each node i: 



2 Otherwise there can be flux through holes or handles of the man- 
ifold. In the graphical expansion, the sum over the resulting flux 
sectors is responsible for killing loop configurations which do not 
correspond to percolation configurations. 

3 However, at this special value of q we can define fc-leg operators 
for odd k by taking O k = S' 1 . . . S' k . 



Z = Tr exp 



E 



Ei 



(24) 



nodes % 



In order to define Ei, we denote the links surrounding i 
by 1, 2, 3, 4, with the weight p pairing being 1 with 3 and 
2 with 4, and the weight (1 — p)q pairing being 1 with 2 
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and 3 with 4: 

exp{-Ei) = p(S 1 .S 3 )(S 2 .S 4 ) + (1 -p)q(S 1 .S 2 )(S 3 .S 4 ) 
+ (l-p)(l-q)(S 1 .S 4 )(S 2 .S 3 ). 

These three terms are in correspondence with the three 
possibilities in Fig. [2] and a graphical expansion of the 
above partition function gives the sum over loop config- 
urations defining the CPLC (Eq. [2]). 

The above Boltzmann weight again has a J. 2 gauge 
symmetry: on changing the sign of the spin on the link 
ij, both e~ Ei and e~ Ej change sign but the overall Boltz- 
mann weight is unchanged. The naive continuum limit 
is again the RP n_1 sigma model described above, and 
the universal properties are expected to be identical with 
those of the IPLC. In Sec. |III G| we will argue that the 
effective fugacity for J. 2 vortices changes sign on the line 
q = 1/2, just as for the IPLC. 



The lattice field theory (24) and the continuum limit 



(21 1 require a couple of comments. Firstly, the PP" de- 



scription does not apply on the boundaries of the phase 
diagram: here the model has a higher symmetry which is 
not taken into account in (24) — see next section. Sec- 



ondly, note that exp(— Ei) is not always positive. This is 
not an obstacle to taking the continuum limit — in fact it 
is an important feature of the problem (note that lattice 
actions which are not real are not necessarily patholog- 
ical, and are the norm in quantum problems). We have 
seen that different signs for the vortex fugacity distin- 
guish the two short-loop phases in the IPLC, and we will 
see that the same is true for the CPLC. (The possibility 
of non-positive Boltzmann weights can important even in 
the Goldstone phase - Sec. IV B). 

In the model without crossings, discussed below, it 
is again crucial that the lattice Boltzmann weight is 
not positive since the continuum Lagrangian contains an 
imaginary 0-term [3] . We now briefly review this formu- 
lation and discuss its implications for the phase diagram 
of the CPLC. 



F. Phase diagram boundaries and CP™ 



As mentioned in Sec. |II A[ the CPLC has an additional 
symmetry [9[ [TT] on the boundaries of the phase diagram 
(when p = 0, or q — 0, or q = 1). On each of the three 
boundaries of the phase diagram, the links of the lattice 
can be assigned fixed orientations, with two incoming and 
two outgoing links at each node, such that the allowed 
pairings are always between an incoming and an outgoing 
link. This orients all the loops. The necessary choice of 
link orientations differs for each of the components of the 
boundary, as mentioned in Sec. |II A| 

Then instead of constructing a lattice field theory us- 
ing real spins S, we may take complex spins z. The inner 
products Si. Si' in the Boltzmann weight for a given node 
are replaced with z} -Zi', where I is the outgoing link and 



V the incoming link. One way of understanding the ap- 
pearance of complex fields is to view the loops as world- 
lines of quantum particles in 1 + 1 dimensional spacetime: 
the fact that we are now dealing with oriented worldlincs 
means that these particles carry a U(l) charge. 

The SO(n) global and J. 2 gauge symmetry in the in- 
terior of the CPLC phase diagram are promoted to an 
SU(n) global and U(\) gauge symmetry on its boundary. 
The appropriate sigma model involves a field on complex 
projective space with a 9 term [SHUT]. The ne ld Q is now 
a traceless Hermitian matrix (Q ab = z a z* b — S ab /n): 



C 



CP" 



= f tr(VQ) 2 



— e MV tr QV M QV„Q. (25) 

Z7T 



The CP™ - field can sustain skyrmion textures, since 
^(CP"" 1 ) = Z. The 9 term weights skyrmions by a 
factor e lB and antiskyrmions by e~ lB . 

Bulk properties of the sigma model depend on 9 only 
modulo 2ir. For n < 2, there is a critical point at 9 = it 
mod27r; other values of 9 are massive (flowing to 9 — 
0mod27r, K = under RG). The 9 = ir critical point is 
the critical point of the loop model at p = 0, q = 1/2 
[21 [TT] . Close to this point the bare value of 9 varies as 



tt) cx (q - 1/2). 



(26) 



For our present purposes the CP™" 1 description tells 
us two things. Firstly, it implies that the left- and right- 
hand boundaries of the phase diagram Fig. [3j which cor- 
respond to the Manhattan lattice loop model, are local- 
ized for all p as previously expected [27J [36] . Here 9 can 
be shown to be equal to zero mod 2ir, so that the sigma 
model is in the disordered phase. However the correlation 
length £, and the typical loop size, diverge exponentially 
as p — > 1, 



e~(i-p) 



-2 g const./(l- 



■P) 



(when q = 0, 1) 



(This follows from the beta function for the CP"" 1 model 
[44] and the fact that the bare stiffness is of order (1 — 
It is the behaviour of £ for non-critical localisation 
in class C [31] , to which the Manhattan lattice loop model 
is related [27].) 

Secondly, the CP n_1 description of the model without 
crossings gives a way of seeing that the vortex fugacity in 
the PP™ -1 description of the CPLC chan ges sign on the 
central line q = 1/2, just as it does at the central point 
q = 1/2 in the IPLC. We now discuss this issue. 



G. Vortex fugacity and the 6>-term for CP' 1-1 

A natural way to approach the field theory for the 
CPLC, at least when the weight p for crossings is small, 
is by perturbing the field theory for the model without 
crossings. In the CP™ -1 language, making the weight p 
for crossings nonzero corresponds to adding a small mass 
for the imaginary part of Q, i.e. 8C oc — ptr (ImQ) 2 |12j . 
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This favours real Q and leads in the infrared to an order 
parameter living on RP™" 1 . 

If we simply set Q to be real, the kinetic term for 
£ C p re -i becomes that of C KP n~i , and the 9 term vanishes. 



Thus we regain the sigma model Lagrangian (21 ) for the 
loop model with crossings. 

However, this does not mean the 9 term plays no role: 
it vanishes only if we neglect vortices. To retain con- 
tinuity of Q in the presence of a vortex, we must al- 
low Q to become complex in the vortex core, and there 
is then a contribution from the 9 term there. (Similar 
mechanisms pertain to quantum magnets described by 
anisotropic sigma models in both 1 + 1 and 2 + 1 dimen- 
sions [45].) 

The vortex core corresponds either to a half-skyrmion 
or to an anti-half-skyrmion, depending on the sign of the 
imaginary components of Q. (This can easily be checked 
for n — 2, when the perturbed CP 1 model is the 0(3) 
model with easy-plane anisotropy.) The effective vortex 
fugacity V is obtained by summing over both possibili- 
ties, leading to 

Voce ie / 2 + e~ ie / 2 . 

Since at p — we have (9 — it) ~ (q — 1/2), the vortex 
fugacity in the RP™ _1 description of the CPLC changes 
sign at q — 1/2, just as the vortex fugacity changes sign 
at q = in the IPLC. 




10 10 2 10 3 10 4 10 5 10 6 



FIG. 8. Two- and four-leg watermelon correlation functions 
Gi and G4 in the Goldstone phase. The fits are to the form 
G fe = C k {lnL/r k f k (k = 2,4), see text. 



initial UV cutoff is A, so that S involves modes with 
wavenumber |fc| < A, then the modes in S are limited to 
|fc| < A for the new cutoff A < A, and cj) contains modes 
in the momentum shell |fc| £ (A, A). Integrating <f> out, 
and working to leading order in K~ x and 7, we obtain 
the RG equations 



IV. THE GOLDSTONE PHASE 



dK 

d7 



2-71 

2tt ' 



d", 
d7 



(30) 



The Goldstone phase shows subtle universal behaviour, 
different from that seen in loop models without crossings, 
which can however be understood in detail. Within this 
phase we can work with the 0(n) sigma model [10] . as 
discussed in Sec. IIIID1 



For most purposes the Lagrangian (22) will be sum 



cient, but to calculate the length distribution in Sec. IV C 
we will need to add a small perturbation, 7, which breaks 
the symmetry from 0(n) to 0(n — 1) x Z2. Writing 



S = (S\S X ), 
where S± is an (n — l)-component vector, we take 



S 2 = 1. 



(27) 



(28) 



We briefly recall one-loop RG results for this model [41)1 
147] . which may be obtained easily using the background 
field method. S(x) is decomposed into a slowly- varying 
field S(x) and rapidly varying fluctuations 4>{x): 



n-l 



S*(x) = ~S\x)^l ~ ^{xf + J2 M^K(x). (29) 

0=1 

The e a (x), a — 1, . . . , n— 1, are a set of vectors orthogonal 
to S(x) (there is a gauge freedom in this choice). If the 



Here r is the RG time: after time r, the new cutoff is 
e~ r A. Again, the important point is that K flows to 
large values in the infrared if n < 2 [10] . 

Note that in two dimensions the higher-order 
anisotropies are as relevant as Sj_ at tree level, but less 
relevant at one-loop order. They will be important for 
the polymer phase diagram discussed in Sec. |VII[ 

We now calculate a range of observables both analyti- 
cally and numerically; details of the numerical procedure 
are given in Sec. |VI[ As we discuss in Sec. |VII[ incom- 
patible hypotheses about the universal behaviour at the 
point p = 1/3, g = l/2 have previously been put forward, 
which is one reason for making a careful comparison of 
numerics and theory in the Goldstone phase. 



A. Correlation functions 

The watermelon correlation function G k (r) may be 
expressed as the two- point function of the operator 
S 1 ... S k (x) (Sec. |IIIB[). Including the UV cutoff A and 



the sigma model stiffness K explicitly in the argument of 
Gfe, a simple calculation following [55] gives 



G k (rA,K) 



(1 + f^lnArr 



(31) 
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where the exponent in the denominator depends on k and 
on the loop fugacity n: 



Oik 



k(k + n-2) 



(32) 



The correlation function in the numerator of (31 ), which 



is evaluated at a separation of the order of the new UV 
cutoff, tends to a constant for large r. Thus the asymp- 
totic behaviour of the watermelon correlation functions 
Gfc(r) is given by a universal power of logr: 



G k (r) 



C k 



(lnr/r ) c 



ro 



A" 



(33) 



with nonuniversal C k , r . 

It is interesting to note that although the stiffness K 
flows to infinity in the Goldstone phase - which we would 
usually think of as implying long range order — all the 
correlation functions G21 decay to zero at large distances 
for n > 0. Correlators G k with odd k have no meaning in 
the RP" _1 sigma model or in the CPLC, but can be de- 
fined in loop models described by the 0(n) sigma model 
without a Z2 gauge symmetry; such models allow oper- 
ators which insert dangling ends. For n — 1, G\ tends 
to a constant, indicating that the entropic force between 
the ends of an open strand inserted into such a soup of 
closed loops vanishes at large separations. 

We may contrast this with the case n = 0, which de- 
scribes the universality class of the dense polymer with 
crossings [TU] — a single loop whose length is compara- 
ble with the number of lattice sites. Here, G2{x) is a 
constant at large separations, in consequence of this fact. 
On the other hand G\{x) has a negative exponent ai, 
indicating that if the polymer is an open strand the two 
ends suffer a weak entropic repulsion. 



We might have expected the logarithmic form of (33) 



to prevent us from seeing the universal exponents a k nu- 
merically, but this is not the case. Fig. H shows G 2 (L/2) 
and G±{L/2) for L x L systems with periodic boundary 
conditions, with L rang ing up to L = 10 6 for G2 and 
L = 10 4 for G4 (see Sec. VI for further details). Simula- 
tions are at p = q = 1/2. We fit G2 and G4 to the form 
Gk = Ck{lnL/r k )~ ak , leading to exponents consistent 
with p2l): 



a 2 = 1.9(1) 



a 4 = 12.5(10). 



(34) 



We have lnr 2 = -15.4(14), lnr 4 = -18(2), consistent 
with the fact that ro is shared between different G k in 
Eq.[33) 
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FIG. 9. The logarithmic increase of the mean spanning num- 
ber (n a ) with system size in the Goldstone phase. Main panel: 
p — 1/2, q — 1/2; lower inset: p — 1/3, q — 1/2. Fits are 
to ^(lnL/Lo + lnlnL/Lo) with lnL ^ -13.78, -8.06 for 
p = 1/2 and p — 1/3 respectively. Upper inset shows the nu- 
merical estimates for the slopes d{n s ) /dlnL plotted against 
InL - they are expected to converge to 1/2-7T (horizontal line). 



This is the number n s of curves which traverse the cylin- 
der from one boundary to the other. Note that n s must 
be even if L is even and odd if L is odd. 

To calculate n s , the correspondence of Sec. (TTlE] be- 
tween the loop model and a spin model must be extended 
to the case with dangling boundary links. We simply take 
the spins on the dangling links to be fixed, with 



■-Hop 



(cos#, sin#, 0, . . . , 0), 



She 



(1,0,. ..,0) 



on the top and bottom boundaries (above we temporarily 
revert to the lattice normalization of S). The graphical 
expansion then goes through as before, except that span- 
ning curves acquire an additional weight cos 9^ 

Denoting the partition function with the above bound- 
ary conditions by Z(9), we therefore have 



((cos 5)"-) 



Z(6) 

z(oy 



(35) 



In the Goldstone phase the stiffness flows to large val- 
ues in the infra-red, so to calculate the right hand side 
we need consider only classical solutions with the appro- 
priate boundary conditions. Letting x be the coordinate 



B. Spanning number 

The logarithmic RG flow of the sigma model stiffness 
K can be seen empirically: this stiffness is directly related 
to the mean 'spanning number' for an L x L cylinder on 
which curves are allowed to terminate on the boundary. 



4 Loops in the interior and curves with both ends on the same 
boundary both retain fugacity n — for example a curve with 
both ends on the upper boundary can either be of colour a = 1, 
in which case it has weight n cos 2 8, or of colour a = 2, when it 
has weight nsin 2 9, and the sum gives n. Spanning curves on the 
other hand must be of colour a = 1, and have weight ncosd. 
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FIG. 10. Probabilities P„ s of n s spanning curves, as a func- 
tion of (n s ) (data for p = 1/2). Curves are the analytical 
expressions from Eq. |38[ with K — (n s ). 



along the cylinder, these are 



S = ±(cos</>(a;),sin(/)(a;),0, ... ,0), <j){x) = 



x(9 + nm) 
L 



Both odd and even values of m are allowed — the bound- 
ary condition is satisfied only up to a sign — but when 
L and m are both odd the Boltzmann weight acquires 
an additional minus sign, as can be seen from the lattice 
partition function^] 

The action of these solutions is calculated using the 
renormalized stiffness K on scale L, leading to 



((cos^^R'^exp 



K 



(9 - irmf 



(36) 



This logarithmic flow (for n = 1) is seen in Fig. [9] for 
two points in the Goldstone phase. We have fitted the 
data for large sizes to the slightly more accurate form 
(n s ) ~ j- (lnL/Lo + lnlnL/Lo) which comes from in- 
cluding the subleading 0(1/ K) term in the beta function 
for the stiffness (23). In the upper inset to Fig. [9] we plot 
the numerical value of the slope d (n s ) /dlnL, which is 
seen to converge slowly to 1/27T for large L. 

Since all cumulants are proportional to K, their ratios 
are universal numbers which we can compare with data: 



These relations are obeyed to good accuracy — plotting 
the two cumulants against (n s ) for p = 1/2 and various L 
gives straight lines with slopes 0.668(5) and 0.274(18) re- 
spectively (data not shown) . Note that the scaling of the 
cumulants means that for very large (n s ) the probability 
distribution of the spanning number becomes Gaussian 
(away from its tails). 

To extract the probability distribution for small integer 
n s we set cos 9 — e in ( 36 1 , 



2 y 



arccos e ) 



+ (-)V 



arccos e — 7r 



(38) 



For even circumference, this gives for example 

P = 2e-^^'\ P 2 = !!M^lM e - 2 (" s >/8. 

In Fig. [10] the expressions for Pq, . . . , Pg are compared 
with data (at p = q = 1/2 and L in the range 10 2 — 
10 6 ) showing remarkable agreement. There is no free 
parameter in these fits. 



For a given value of 9, only one or two values of m are 
not exponentially small in K . 

To extract low-order cumulants for n s , we set cos 9 = 
e~ x and expand in x. Since (36 1 is dominated by the 



m = term for 9^0, the difference between even and 
odd L is not seen. The Ith cumulant is given by 



K 



2 ^(arccos e*) 2 \ y=0 . 



In particular, the mean spanning number is given by the 
renormalized stiffness K, so 



(n s ) 



2 - n , L 
— — In — . 
2tt L 



(37) 



In the mapping of Sec. |III E| this sign comes from factors of 
Stop-Si and Sbottom-S; in the Boltzmann weight for links I ad- 
jacent to the boundary links. 



C. Length distribution 

To calculate the length distribution for a loop we must 
consider the RG flow away from the Goldstone phase 
induced by a symmetry-breaking perturbation. To sum- 
marize the result of the following calculation, which is for 
n = 1, the probability for a loop randomly chosen from 
the soup to have length I falls off as 



P(l) cx 



1 



Pln"(l/l ) 



(39) 



for large I. 

Fig. |ll| shows the distribution obtained numerically for 
loops of length up to Z ~ 10 10 . We multiply P(l) by I 2 in 
order to expose the logarithmic correction, which we fit 
to the form a(\nl/lo)~ c . We obtain 



c= 2.03(3) 



(40) 



in striking agreement with ( 39 1 
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FIG. 11. Main panel: The probability distribution P(l) for 
the length of a loop in the Goldstone phase. We have mul- 
tiplied by I 2 to remove the expected power law part, leaving 
the logarithmic dependence. The fit is as described in the 
text, with ln/o = —33.7(8). Data is for p = q — 1/2. Inset: 
two ways of defining the effective finite-size exponent value 
TeB (see text). Green squares: data from P(l), together with 
the fit (dashed line) implied by Eq. [39](lnZo = -32.7). Red 
circles: data from AX (I). 



Note that P(l) differs by a factor of I from the length 
distribution for the loop passing through a fixed link, 



P 



fixed 



itov(Z) cxlP(l), 



simply because longer loops visit more links, 
evaluated using P(l) is finite, as for (39). 



(41) 
Thus (I) 



Let g(x) denote the generating function for the length 
of a loop randomly chosen from the soup (the 'marked' 
loop) : 

g{x) = p ( l ) xl = (z lcnsth ° f markcd lo ° p ) • (42) 



In order to extract g(x) we use the trick of Ref. [48] . 
splitting the components of S, or equivalently the loop 
colours, into two groups. For simplicity we will consider 
only the loop model at fugacity one, though it would 



be easy to generalize. We split S as in Eq. |27j yielding 
an (n — l)-component vector S±. In the graphical ex- 
pansion of a lattice model, say of the CPLC (the IPLC 
would be similar) we correspondingly split the loops into 
unmarked loops, whose colour index is equal to one, and 
marked loops, whose colour index runs over 2, . . . , n. Af- 
ter summing over loop colours, a configuration with N 
marked loops has a weight (n — 1) N , and expanding the 
partition function in (n— 1) is equivalent to expanding in 
the number of marked loops in the configuration. Writing 



1. 



Z{n')=J2 W C+n' Wc 



(43) 



C; one 
marked loop 



Here Wc is the weight of a configuration C in the CPLC 
at n = 1 (Eq. [I]). The first term on the right hand side 
is equal to one. 

Next, we wish to modify the weight of a configuration 
by the factor x , where I is the total length of marked 
loops. This is easily done: all the inner prod ucts S t. Si; 
appearing in the node factors e~ Ei (see Sec. HIE) are 
replaced according to: 



Si. Si' 



Si Su + x Si±.Si'±. 



(44) 



In this way every unit length of marked loop picks up a 
factor of x, and the graphical expansion yields 



Z{ri, x) = l + n' 



E 

C; one 
marked loop 



We x 



(45) 



Differentiating with respect to n' gives the required gen- 
crating function: 



d„> Z(n',x)\ n>=0 



(46) 



d n , Z(n',l)| n#=0 
In terms of the free energy density f(n', x), 

Z(n>, x) = e - L ' 2 ^ n '' x \ f(n',x) = f + n' h(x) + ..., 

this is g(x) = /i(x)//i(l). In order to calculate the 
length distribution for large values of we require the 
free energy for x ~ 1 and to first order in n' . 

In the continuum description, the symmetry-breaking 
perturbation x 7^ 1 leads to an infinite number of relevant 
perturbations of which the most relevant is that in Eq. |28| 
Setting x — exp(— fj,) with /1 <^ 1, we have 7 ~ fi, the 
constant of proportionality being nonuniversal. 



Beginning with the Lagrangian of Eq. 28 we integrate 
out high frequency modes, retaining their contribution to 
the free energy, up to an RG time . This gives 



/(if, 7) = f < (K, 1 ) + f > (K, 1 ), 



(47) 



where we have split up the contribution from the modes 
that have been integrated out, 



/<(^,7) 



n 
4tt 



T *^(lnX(r)+ 7 (r))-n'A 



(the nonuniversal constant A ensures / < 0, as required 
by Eq. 45 1 and those remaining: 



f > {K, 1 ) = e- 2T 'f{K{r,), 1 {r,)). 
The solutions to the RG equations |30|) are 



K(t) = K 



2tt' 



7 (t) = 7e : 



iT / 2-kK 
V 2nK + 



(48) 



Stopping the RG when 7* becomes of order one, f > may 
be approximated as the free energy of a massive theory 
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in which S executes only small fluctuations around S 
(1,0,. ..,0): 



C = ^ ((V^) 2 + 7 *5i + 0(Si, 



(49) 



(The 0(S 2 _/ K*) term comes from the sigma model mea- 
sure.) However, the dominant terms come from f < . 
the order we require, 



To 



f{K, 7) = —B + 7 ( 2nK 



8n 2 K 2 
In I/7 



(50) 



(B is a constant.) We thus have the form of the generat- 
ing function at small fi: 



(e-» l ) = l-C f x 1 



AttK 
lnl//i 



(51) 



The constants C and K are non-universal, and the fact 
that the leading \x dependence is linear in 11 is simply 
a consequence of the fact that (I) is finite. However we 
may infer the behaviour of P(l) at large / from the form 
of the nonanalytic term, yielding Eq. |39| 

For a general critical loop ensemble, P(Z) ~ l~ T for 
some t > 2, and the mean size AX of a loop scales with 
its length as AX ~ l x l d s . The fractal dimension dj is 
related to r by 



r = 2/d f + 1. 



(52) 



In the Goldstone phase, t = d e g — 2, with logarith- 
mic corrections. We may define finite size estimates of r 
either using dlnP/dlnZ or using dlnAA/dlnZ and the 
scaling relation. These are plotted in the inset to Fig. 11 
both are expected to converge to two, but with different 
logarithmic corrections. 

Here, AX is defined as the mean extent of a loop in one 
of the coordinate directions. A similar quantity — the 
mean square end to end distance of an open trail in the 
ISAT model — was considered numerically by Owczarek 
and Prellberg |16j . and logarithmic corrections to Brow- 
nian scaling were found. It would be interesting to cal- 
culate these quantities analytically. 



V. THE CRITICAL LINES 

The critical lines separate the Goldstone phase from 
phases with short loops. In the language of the RP™ -1 
model, they correspond to order-disorder transitions at 
which Z2 vortices are set free. Here we give numerically 
determined critical exponents for this transition at n = 
1 and briefly consider an approximate RG treatment of 
vortices [151120]. 




0,86 



FIG. 12. Main panel: the mean spanning number n 3 as a 
function of q for p = 1/2 and various system sizes, showing a 
crossing at the transition. Inset: data collapsed according to 
Eqs. p| [55l 
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FIG. 13. Main panel: the behaviour of the mean spanning 
number at p = 1/2 very close to the critical point. Note the 
larger system sizes compared to Fig. |12| Inset: finite size 
corrections to the spanning number at the critical point (note 
log- log scale) and linear fit leading to estimate of j/i rr . 



This is manifested in the crossings of the various curves 
in Fig. [12] which shows the spanning number Tig as a 
function of q for p = 1/2 and for cylinders of various 
sizes. 



Fig. 13 shows the same quantity, but in the im- 
mediate vicinity of the critical point and including much 
larger system sizes (up to L — 128, 000). The main panel 

0.3; here finite size effects 
are much stronger. 



of Fig. 14 shows data for p 



- visible in the drift in crossings 
The basic finite size scaling form for n s is 



A. Critical spanning number, u, and y\ r 



n s = h(x), 



(53) 



At the critical point, the dimensionless quantity n s (de- 
fined in Sec. IV B I is expected to take a universal value. 



where 5q = q — q c . We take into account also nonlinear 
dependence of the scaling variable x on dq 7 replacing the 
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FIG. 14. Main panel: spanning number at p = 0.3 (close-up 
of the critical point) showing larger drift in crossing points 
than at p = 1/2. Inset: the vertical coordinate n* is the 
crossing in the spanning number between consecutive system 
sizes and L* is the geometrical mean of the sizes. Red/black 
dots correspond to p = 1/2/p = 0.3. The horizontal line is 
the estimated asymptotic value ( 56 1 . 



FIG. 15. Main panel: the two-leg watermelon correlator Gi 
at the critical point. The power-law decay Gi ~ L~ 2x2 gives 
the fractal dimension df via ( |60[ ). The upper inset compares 
this value (indicated by the horizontal black line) with the 
finite-size estimates coming from AX (I) (red circles) and P(l) 
(green squares) — see text. The lower inset shows the four-leg 
watermelon correlator G4. 



second equation above with 

x = L 1/v 5q (1 + p x 8q + I3 2 5q 2 ) 



(54) 



and finite size corrections with (negative) irrelevant ex- 
ponent 2/i rr in the form: 



n s = h{x){l + Ly-{k+p A x)) 



(55) 



A reasonable scaling collapse may be obtained by ad- 
justing the values of q c , Pi, v and y xrc . To find these 
values we fit n s to the form (55), constructing h(x) us- 
ing B-splines with 22 degrees of freedom. The result for 
p = 1/2 is shown in the inset to Fig. [l2] What is plotted 
is nf = n s j (1 + L Vi "((3 3 + P^x)), which should be equal 
to the scaling function h(x) by | |55| ). 

Our estimates of the correlation length exponent and 
(universal) critical spanning number, obtained from the 
data at p = 1/2, are: 



2.745(19) 



2.035(10). (56) 



We cannot constrain the irrelevant exponent very pre- 
cisely. From the full fit, we obtain 



Wn 



-(0.2,0.35). 



(57) 



A direct estimate from the finite size corrections to the 
spanning number at the critical point gives a result com- 
patible with this: the fit in the inset to Fig. 13 corre- 
sponds to J/i rr = —0.272. 

Results for p — 0.3 are consistent with our expecta- 
tion that all points on the critical line are in the same 
universality class, but error bars are larger because of 



larger finite size effects and smaller system sizes. We 
find v = 2.87(10) and nf* = 2.07(3). With regard to 
the convergence to a common nf*, see also the inset to 
Fig. [14] which shows the vertical coordinates of the cross- 
ings between curves for consecutive L values. 



B. Watermelon correlators, fractal dimension, and 
length distribution 

Next we consider the watermelon correlation functions 
G2 and G4 defined in Sec. |IIIB[ We evaluate these cor- 
relators at separation L/2 for a range of system sizes L - 
see Fig. 15 The data is for the critical point at p = 1/2. 



Fitting to pure power laws, 

G fc (L/2)ocL- 2 ^ 



(58) 



we obtain the scaling dimensions of the two- and four-leg 
operators: 

x 2 = 0.091(1), x 4 = 0.491(1). (59) 

The scaling relation 

d f = 2-x 2 (60) 

gives the fractal dimension of the critical loops, 

d f = 1.909(1). (61) 

We may obtain independent estimates of df from AX (I), 
the mean linear size of a loop of length /, and from the 
length distribution P(l) and the scaling relation (p32|. 
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The finite size estimates for df coming from the numerical 
estimates of din AX/dlnJ and dlnP(Z)/dlnZ are shown 
in the inset to Fig. [15] Both plots are consistent with 
(611. 



yields at n = 1 



C. RG equations in the presence of vortices 



For an approximate description of the transition, we 



extend the RG description (23) to take account of the 
nonzero fugacity for Z2 vortices. In this we follow the 
treatment by Fu and Kane of the 0{2N)/0(N) x O(N) 
sigma model at N — > 20 ]. This sigma model and the 
RP™ _1 sigma model are similar — both sustain Z2 vor- 
tices, and each reduces to the XY model in an appropri- 
ate limit, which can be expanded around. An expansion 
around the XY limit was also considered for the 0(n) 
model near n = 2 in Ref. |49) . The importance of topo- 
logical defects in replica sigma models for localisation in 
two dimensions was also pointed out by Konig et al. , who 
developed an RG approach to localisation in the chiral 
symmetry classes taking account of Z vortices [19] . 

Since RP 1 = S 1 , our sigma model coincides with the 
XY model at n — 2. With the normalisation of Eq. [21] 
this has a Kosterlitz Thouless transition at the critical 
stiffness K c = 8/77. The RG flow near this point is gov- 
erned by the Kosterlitz RG equations for K and the vor- 
tex fugacity, which we denote V. 

We assume that we can expand the RG equations in 
(2 — n), and that V should be interpreted as the fugacity 
for J.2 vortice^ when n ^ 2. At lowest order, these 
equations are corrected by the /3-function for K in the 
absence of vortices, dK/dr ~ (2 — n)f(K): 



dV 

"d7 

dK 
dr 



= (2 - n)f(K c ) - V 2 



(62) 
(63) 



These equations yield critical points at K = K c = 8/tt 
and V = ±-\/ (2 — n)f(K c ), with critical exponents 



(2 - n)Trf(K c 



Vh 



(2 - n)wf(K c ) 



The critical stiffness K c ~ 2.4 is of roughly the same 
magnitude as the critical winding number n c s nt ~ 2.0 
at n = 1 (56), as we expect from Sec. IV B Making 
the further approximation of evaluating f(K c ) using the 
perturbative /3-function at large K, 



f{K) 



2?r 



2ttK 



For n > 2, a Z2 vortex corresponds to e.g. S = 
(cos 9/2, sin 9/2, 0, 0); note that when n = 2 this configration 
becomes a Z vortex of charge ±1. (Here 9 is a polar coordinate; 
recall S~ -§.) 



1.9, 



2/irr 



-0.5. 



(As n — > 2, v diverges and the irrelevant exponent tends 
to zero.) 

As expected, this crude approximation does not give 
quantitatively accurate results for n = 1, but it does 
reproduce the qualitative structure of the phase diagram, 
with the Goldstone phase sandwiched between massive 
phases at positive and negative V, and the appearance 
of a large correlation length exponent. 

A comparison with alternative approaches — for exam- 
ple an approximate treatment of sigma model directly at 
n = 1 in the supersymmetric formulation, or an expan- 
sion in (2 — n) avoiding the additional large K approxi- 
mation required here — would be desirable. 



VI. NUMERICAL METHODS 

For system sizes up to L ~ 2 x 10 4 we can use a straight- 
forward Monte Carlo procedure, which of course benefits 
from the fact that node configurations are independent 
random variables when n = 1. Very large sizes require a 
more efficient 'knitting' procedure. 

In the straightforward approach, we construct indepen- 
dent L x L samples and calculate the spanning number 
n s , the length distribution P(l), the average loop size for 
a given length AX (I), and the correlation functions G2 
and G4 at separation L/2. G2 is the probability that 
two links lie on the same loop, and G4 is taken to be the 
probability that two nodes are visited by the same pair 
of loops. This is not the only way that two nodes can be 
connected by four strands - the four strands can be part 
of a single loop - but the scaling is the same for both 
types of contribution. 

The data for P(l) and AX (I) are stored in histograms 
in a logarithmic scale in I; to avoid finite size effects we 
discard boxes whose AX is within two standard devia- 
tions of L from fits. 

Note that when n = 1 the Boltzmann weight does not 
depend on the boundary conditions, so we can go trivially 
between the fully periodic BCs used to define Gk and 
the cylindrical BCs used to define n s . We have also used 
another geometry, a square with reflecting BCs for the 
loops on the edges but links sticking out at the corners, 
to measure the relative probabilities of the three distinct 
corner connectivities as a function of p, q, and L. This 
gives a simple means of determining the phase diagram 
in Fig. § 



A. Knitting 

As we expect logarithmic behaviour in the Goldstone 
phase, it is crucial to be able to study very large sys- 
tems. The above procedure is very efficient in terms of 
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CPU time, but it is limited by the available computer 
memory since it requires us to store the full configura- 
tion. The first improvement comes from the fact that a 
cylinder of circumference L can be 'knitted' by the suc- 
cessive addition of L x V strips, where L' is of order one. 
At a given stage in the growth of the cylinder, it con- 
tains both closed loops in the interior, and strands which 
end on the boundary. The only information about the 
configuration which we need to store is the connectivity 
of the links on the boundary, together with the lengths 
of the connecting strands. (Note that this connectivity 
information is configuration-dependent — this is distinct 
from the transfer matrix approach, in which the transfer 
matrix is not configuration-dependent, and which is lim- 
ited to small sizes.) When we add a new strip some loops 
will become closed: we add their lengths I and their sizes 
AX - defined as the height of the loop in the growing 
direction - to our histograms. 

Growing the cylinder layer by layer, the required 
memory is proportional to the circumference L (strictly 
L\nL) rather than to L 2 as in the straightforward ap- 
proach (algorithms for percolation which avoid storing 
the full configuration also exist [SU]). We have been able 
to compute cylinders with L up to 10 6 and height much 
greater than L. 

This method requires storing the (realisation- 
dependent) connectivity information for the boundary 
links of a large cylinder. This has the flavour of a trans- 
mission matrix in a localisation problem. It would be 
interesting to consider the properties of this 'matrix' in 
more detail. 



VII. POLYMER COLLAPSE 

A long polymer with repulsive (or excluded volume) 
interactions between segments displays the universal be- 
haviour of the lattice self-avoiding walk (SAW), while 
strong enough attractive interactions cause the polymer 
to collapse. The boundary between these two regimes is 
the so-called point. In de Gennes' description of the 
polymer via the 0(N — >• 0) model, the SAW corresponds 
to the critical point, and the point to the tricritical 
point [U 021 [51] . However the actual situation is more 
complicated, especially in two dimensions. 

While the SAW behaviour is extremely robust, the 
point is more subtle. Different lattice models can yield 
different universality classes of collapse transition [T6l HH 
1521 [S3], and it also turns out that the phase diagram 
in the vicinity of the point can have a more compex 
structure than would naively have been expected [HI 152] 

The interacting self-avoiding trail model (ISAT) has 
been particularly controversial (TBI [HI [2H EH] [57] , with 
numerous conflicting results and hypotheses put forward 
for the critical exponents at the collapse point. In this 
model, the polymer can visit links only once, but nodes 
twice, so allowed configurations of a closed polymer are 
equivalent to allowed configurations of a single loop in 
the CPLC. Here we give a field theoretic description of 
the ISAT which explains the phase diagram found nu- 
merically [18] and shows that the ISAT point is highly 
fine-tuned from the point of view of more general poly- 
mer models, being an infinite-order multicritical point at 



B. Shuffling 

A further improvement that significantly reduces CPU 
effort, and allows us to calculate G2, involves construct- 
ing sets of strips of width L and height H = L/20. 
For each strip we store the connectivity of the boundary 
links. Joining 20 of them yields a square sample, and we 
have enough information about the sample to construct 
G2(x,y) (for links x and y which lie on the boundaries 
of two strips) as well as n s . From each set of 20 pieces, 
many different samples may be created by shuffling the 
order of the pieces and by rotating them in the trans- 
verse direction. We construct 1000 samples for each set. 
These samples are not of course independent, so we es- 
timate error bars by producing many independent such 
sets (80 — 200) and examining the statistical flucutations 
between sets. 

For the calculation of n s in large systems to a given pre- 
cision, the shuffling and rotation of strips reduces CPU 
time by a factor of 200. 

We have parallelized this procedure for OpenMP and 
for CUDA GPUs. The parallelization for the graphics 
cards was particularly efficient - a typical programme 
ran almost 100 times faster on an Nvidia Tesla M2070 
card than on a single core in an Intel Xeon E5520 CPU. 
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Crystalline phase 
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Liquid phase \ 
<S^>*0, <S'>*0 



Zero density phase 
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t (interaction parameter) 

FIG. 16. Schematic representation of the phase diagram for 
the polymer found in Ref. [18] . together with our interpreta- 
tion in terms of a perturbed 0(1 + n') model. The thin solid 
line is in the universality class of the self-avoiding walk (poly- 
mer in good solvent). The dot is the multicritical collapse 
point, or Q point, with full 0(1 + n') symmetry. The bold 
line is a first-order transition. The dashed line is a transition 
in the Ising universality class between two regimes in which 
the polymer is dense (visits a finite fraction of the links of the 
lattice) - the 'crystalline' and 'liquid' phases. 
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which the symmetry of the problem is enhanced. 

On a lattice of large finite size L, the ISAT partition 
function iifl 



E 

polymer 
configs 



^lcngth^no. self-contacts 



(64) 



A self-contact is a node visited twice by the polymer. We 
are now at fixed length fugacity rather than fixed length, 
but the two ensembles are simply related^] t mediates 
interactions which are repulsive for t < 1 and attractive 
for t > 1. 

An interesting feature of this model is the phase dia- 
gram, obtained numerically by Foster [18] . A schematic 
version is shown in Fig. [l6j For small k the polymer is 
of finite typical size and not critical (the 'zero density' 
phase), while for large k it is dense, i.e. has length of 
order L 2 . When t is small the transition between these 
phases shows the usual critical behaviour of the SAW 
(thin solid line), and for large t there is a first order 
transition associated with the collapsed polymer (thick 
line). The point separates these. 

An unexpected feature, from the point of view of the de 
Gennes theory, is an additional line of transitions within 
the dense phase (dashed line), which is found numerically 
to be in the Ising class [IS] . A similar line of Ising tran- 
sitions terminating at a 9 point was found in a polymer 
model without crossings studied by Blote and Nienhuis 
|52j , and a heuristic explanation was provided by associ- 
ating Ising degrees of freedom with the faces of the lattice 
[52] [58] , for which the polymer was a domain wall. 

For a field-theoretic description of the ISAT, we make 
use of the fact that precisely at the 8 point, which is at 
k = 1/3, t = 3 QI], the ISAT maps to the CPLC at a 
point in the Goldstone phase, namely n — 1, p = 1/3, 
q = 1/2. The point is therefore described by the 
0(n —> 1) sigma model [TO], and we can understand the 
region around it by perturbing this sigma model. The 
following considerations may readily be generalised to 
the various modifications of ISAT that have been studied 
numerically, e.g. on other lattices [51], with additional 
interactions [551461] . or in higher dimensions |6"2"] . 

Before continuing, we note that an alternative conjec- 
ture for the critical behaviour of the ISAT point was 
put forward on the basis of numerical transfer matrix 



7 Translated versions of a loop are counted as distinct polymer 
configurations in the above sum. 

8 Fixed length fugacity k is more convenient from the point of 
field theory and also if we wish to consider phases in which the 
polymer becomes space-filling. The two ensembles are related by 
a Laplace transform. The critical properties of a polymer of fixed 
large length on an infinite lattice are governed by the singularity 
of the partition function closest to k = 1421 — i.e. in the case of 
Fig. |16| by a point on the line separating the zero-density phase 
from the dense phase. The 'crystalline' and 'liquid' phases are 
dense, i.e. the length of the polymer is comparable with the 
lattice area. 



calculations in Ref. [15] ■ According to this conjecture, 
the point has nontrivial critical exponents identical to 
those of an exactly solvable model of polymers without 
intersections [521 ESS] ■ This is at odds with the predictions 
of the Goldstone phase, which yields trivial critical ex- 
ponents together with universal logarithmic corrections. 
We believe the Goldstone phase scenario is convincingly 
established by our numerical results for large systems 
(note that Sec. IV B includes data at the relevant point 
p = 1/3, q = 1/2), together with the logarithmic be- 
haviour seen numerically in Refs. [TBI fT7] and the theo- 
retical arguments of Ref. [TU] and Sec. |III[ and that the 
apparent nontrivial exponents in Ref. [18] are due to log- 
arithmic finite size correct ions 

To describe the single-polymer problem in the lan- 
guage of the spin model, we proceed along similar lines to 



Sec. IV C expanding the partition function for the CPLC 
in n = n — 1 to separate out configurations with a single 
marked loop (which will be our polymer). 

To control k and t for this polymer we must modify the 
partition function ( 24 ) in a way that breaks the symme- 
try of the 0{n) sigma model down to Z 2 x 0(n'). At each 
node, (24) contains terms of the form {S\.S2){Sz-Si). 
Such a term corresponds to two sections of loop pass- 
ing through the node, one connecting link 1 to link 2, 
and one connecting link 3 to link 4. In the new ensem- 
ble, these sections can be sections of marked or unmarked 
loop, and we modify the weights accordingly: 

(S ? 1 .^ 2 )(S ? 3 .^ 4 ) — ► (66) 
{S\Sl){SlS\) 
+ 3k(Sl.S^)(S ±3 .S ±4 ) + 3k(S ±1 .S ±2 )(SlSl) 
+ 3k 2 t(S ±1 .S ±2 )(S ±3 .S ±i ). 

Each unit of marked length acquires a factor of 3fc, and 
each meeting of two marked strands acquires an addi- 



tional factor of t/3. Expanding in n' as in Eqs. 43 45 
Z(n',k,t) = 
1 



n 



E 

C; one 
marked loop 



W c (3fc) length (t/3) r 



self-contacts 



+ 



9 We can use the sigma model to estimate the eigenvalues of 
the transfer matrix for a cylinder of circumference L. Coarse- 
graining by a factor L fields an effective ID sigma model with 
renormalised stiffness K^. This is equivalent to the quantum 
mechanics of a particle on an (n — l)-dimcnsional sphere, with 
energy levels 

_, _ fc(fc-l) 



0,1,2, 



(65) 



2K L 

in the limit n — > 1. This leads to the finite size estimates 
xj. = Efe/27r for the scaling dimensions, which tend to zero log- 
arithmically with L. For a very crude estimate at small sizes, 
we can approximate by the spanning number on an L X L 
cylinder, which we have calculated numerically at p = 1/3. For 
L from 5 to 12, this gives X2 varying from 0.103 to 0.086. 
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Separating the sum into a sum over configurations of the 
marked loop (polymer) and a sum over the configurations 
of the other loops, and performing the latter, 



Z{n',k,t) = 1 + n' Z po} (k,t) + .. 



We discuss only the partition function, but one may eas- 
ily check that the natural geometrical correlation func- 
tions in the polymer problem, i.e. the watermelon corre- 
lators, can be expressed as correlators of S± in the replica 
limit n' — > 0. 

The polymer multicritical point at (k, t) — (A, 3) corre- 
sponds to the CPLC at p = 1/3, q = 1/2 and therefore to 
the sigma model with full 0(n) symmetry. Varying (k, i) 
away from this point introduces symmetry-breaking per- 
turbations, of which the most relevant are 71 and 72: 



C = 



K 



((V5) 2 



(68) 



o ±1 = S_ 



"2 



ox, = (sir - 2 -^si 



Before considering the more detailed RG picture, we 



identify the phases in Fig. 16 with the phases of this 
perturbed sigma model. 

These are characterized by whether the Z 2 and 0(n') 
symmetries are broken, i.e. whether (S 1 ) 7^ and 
whether (Sx) ^ 0. When ^ 0, the polymer fills 

the system densely and the transverse modes are in a 
Goldstone phase. To see that such a phase is possible, 
note that when fluctuations in S 1 are massive we may 
imagine integrating them out to get an effective 0(n') 
sigma model for S±. In the limit n' — > 0, this sigma 
model has a Goldstone phase as a consequence of the 



beta function in Eq. 23 This is just the theory for the 
dense polymer with crossings studied in Ref. 10J. 

The phase with (S 1 ) = (S±) = does not appear 
upon perturbing around the point, since the latter is 
controlled by the fixed point at infinite stiffness. However 
the other three do: 

Zero-density phase, (S±) = 0, (S 1 ) ^ 0. The lead- 
ing effect of reducing k from the O-point value is to intro- 
duce a mass for S±, so that S orders in the longitudinal 
direction and the transverse modes S± are massive. Cor- 
relators thus decay exponentially for the polymer, and it 
has a finite typical size. 

Dense phase with Ising disorder, (Sx) 7^ &: 
(S 1 ) = 0. Moving away from the 0-point by increasing 
k makes S 1 massive, and S 'orders' in the transverse 
plane. The system thus goes from the Goldstone phase 
of the 0(n — > 1) model to the Goldstone phase of the 
0(n — > 0) model. The two- leg watermelon correlation 
function is a constant at long distances, as discussed in 
Sec. |IV A| meaning that the polymer is dense. 

Dense phase with Ising order, (Sx) 7^ & (S 1 ) 7^ 
0. For appropriate values of the coefficients 71, 72, the 
(renormalised) free energy is minimised when both (Sx) 



and (S* 1 ) are nonzero, so that both symmetries, Z 2 and 
O(n'), are brokenp*] 

The field theory description also determines the nature 
of the phase transitions. First, consider the thin solid 
The field S 1 is Ising-ordered on both 



(67) line in Fig. 16 



sides of this transition, and its massive fluctuations play 
no role in the critical behaviour of S±. We therefore 
have the critical point of the 0(n') model in the limit 
n' — > 0. This is the usual description of the self-avoiding 
walk [2 142) , confirming what we expect and find in the 
polymer problem. 

Next, consider the dashed line in Fig. 16 Here, S 1 
undergoes an ordering transition at which Z 2 symmetry 
is broken. Thus we would expect an Ising transition. We 
must check, however, that the massless degrees of free- 
dom associated with Sj_ — which are in the Goldstone 
phase — do not modify the Ising critical behaviour. But 
it is easy to see they do not. The most relevant coupling 
allowed by Z2 x 0(n') symmetry is via the product of the 
energy operators, 



Rising x Eg, 



oldsto 



(69) 



This composite operator has dimension (length) , so is 
irrelevant. 

Finally, consider the thick line in the figure, which sep- 
arates phases breaking different symmetries (Z2 on one 
side, and 0(n') on the other). According to Landau the- 
ory this transition should be first order, as it is numeri- 
cally found to be [T5] , 

What does Ising order/disorder mean for the polymer? 
Define a new configuration of Ising spins fi p on the faces 
F of the square lattice by the requirement that the poly- 
mer is the (only!) domain wall in this configuration. 
Then the dense phase with Ising disorder corresponds 
to antiferromagnetic order in /z, while fi is disordered 
in the dense phase with Ising order. (To show this, we 
write (fiFfiF'} in terms a correlator of twist fields in the 
RP"- 1 model, which force the fields Q la with a > 1 to 
change sign on a line connecting F and F' .) Antifer- 
romagnetic order in /1 is equivalent to the 'crystalline' 
order of Ref. [H] , and it becomes perfect when the poly- 
mer visits every link of the lattice. It is also essentially 



10 Note that the labelling of the phases above is in terms of the ex- 
pectation values/broken symmetries of the perturbed 0(n) sigma 
model, not of the perturbed loop model/RP n_1 model. These 
are not the same, since only gauge-invariant expectation val- 
ues are meaningful in the latter. In particular, the zero-density 
phase does not beak any of the symmetries of the perturbed 
R pn-l mo del. These are Z2 and SO(n') — Q 11 is invariant 
under both, Q la with a > 1 is invariant under neither, and Q ab 
with a, b > 1 is invariant under Z2 only. No expectation value 
breaks Z2 without also breaking SO(n'). (In the zero density 
phase the unmarked loops of course remain in the Goldstone 
phase of the 0(n — > 1) model — in order to write down the 
associated watermelon correlation functions, we would have to 
extend S 1 to a vector Si with m — > 1 components, increasing 
the symmetries of the model.) 
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equivalent to the Ising order denned in Refs. [52 , [SB] for 
a different model. 

For a more detailed picture of the phase diagram, we 
use the RG equations for the sigma model. We must 
include the lowest two anisotropics as in Eq. 68 where 
7i j 72 are linearly related to the perturbations 



Sk = k- 1/3, 



5t = t-3 



(70) 



when these are small (we will give approximate expres- 
sions below). After running the RG up to a large time 
t*, we have (for n = 1) 



A 



2*. 
2^ : 



7i* 



7i e 



2r» 



(tJ2ttK) 2 



72* 



72 e 



For generic small initial values (71 , 72 ) , we will renor- 
malise to a regime where 71 = O(l) and I72I -C (71 1, 
putting us deep within one of the phases - either the 
zero density phase or the dense phase with Ising disor- 
der, depending on the sign of 71 . The phase transitions 
occur instead in the regime where the renormalised 71 
and 72 become of order one simultaneously. 

Since the stiffness of the renormalised sigma model is 
large (r* ~ ln|7i| -1 / 2 ~ In |t*j| -1 ' 2 ) we may determine 
which phase it is in simply by minimising the potential in 
the Lagrangian, i.e. by mean field theory. In doing this 
we must bear in mind the constraint S 2 — 1. We find the 
three phases described above, with the phase transition 
lines located at: 



SAW : 71 
Ising : 71 
1 st order : 71 



5 72 


( 4ttK 


Vlnl/M 




( AitK 


V lnl /l72| 


b 2 


( AttK 





(72 > 0) 
(72 > 0) 
(72 < 0). 



(These formulas are valid asymptotically close to the 
point.) For a very crude estimate of the relation between 
(71 , 72) and (5k, St) we can evaluate the right hand side of 
Eq. 



GG 



for spatially constant Si, and take the logarithm of 
the Boltzmann weight for a node to obtain the potential 
terms in the Lagrangian. We find 



7i 



C[Sk+ —St 
45 



"2 



-C 



St 
18' 



(71) 



where C is an undetermined constant. 

The above confirms that the three phases meet at the 
0-point, and shows that the SAW and Ising critical lines 
are asymptotically parallel as they approach the point. 

A remarkable consequence of the above field theory 
mapping is that the point of the ISAT - despite the 
simplicity and naturalness of this model - is in fact an 
infinite order multicritical point! We have mentioned the 
two most relevant anisotropies 0±i and 0±2, but there 



is an infinite number of these, with 0±k = {Sj_) k + 
At n — 1, the RG equation for 7^ is 



dTfc = / _ 2fc 2 - k + 1 
dr V 2ttK 



7fc 



(72) 



so all of these are relevant at the K = 00 fixed point. 

The point of the ISAT does not therefore repre- 
sent the universality class of the generic point polymer 
model with crossings. (This explains why the model in 
Ref. [M] shows different behaviour to the ISAT.) A simi- 
lar argument explains why the three-dimensional ISAT 
|62j shows distinct universal behaviour from standard 
models of polymer collapse. 

In two dimensions the generic 0-point behaviour in 
the presence of crossings is different from that [53] in 
their absence. Since the ISAT point does not repre- 
sent the generic behaviour of the 0-point polymer with 
crossings (which we expect to be described by the tri- 
critical 0(n — > 0) model), exact exponents for the latter 
are still unknown. We will discuss RG flows for 0-point 
polymers in detail in a separate publication. 



VIII. OUTLOOK 

The loop models we have discussed are described by 
'replica' sigma models of the kind familiar from localiza- 
tion and polymer physics. Critical points in such prob- 
lems correspond to logarithmic conformal field theories 
[3 3) , which in general have not yet yielded to exact treat- 
ments. We hope that the transitions in the loop models 
will provide a testing ground for new approaches. 

On the analytical side, while exact results for the crit- 
ical behaviour would be desirable, the development of 
more accurate approximate techniques would also be use- 
ful. Numerical worl|^"| on the critical loop model should 
also be extended to other values of n, either via Monte 
Carlo or the transfer matrix 23 , in order to understand 
the whole family of critical points for < n < 2. We plan 
to return to these issues. (Three-dimensional IRP' 1 ^ 1 loop 
models also exist — we will report numerical results else- 
where.) 

The connection between the CPLC at n = 1 and disor- 
dered fermions remains an open question. To begin with, 
recall the situation for the loop model without crossings 
(p = 0). This can be related to localisation in at least 
two ways. Firstly, a limiting case of the Chalker Codding- 
ton model for the quantum Hall effect |65j . in which the 
scattering matrices at a node become 'classical', yields 
the loop model without crossings — i.e. classical perco- 
lation. This is the familiar semiclassical description of 
the quantum Hall transition [66], but because quantum 



11 The IPLC is also a candidate for numerical simulations: we ex- 
pect universal behaviour to coincide with that of the CPLC, but 
this remains to be confirmed. 
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tunnelling has not been taken into account, it does not 
correctly capture the universal critical behaviour. How- 
ever, the loop model has a second relationship with lo- 
calisation which is less obvious and which does not rely 
on suppressing quantum tunnelling. This is due to an ex- 
act mapping from a network model for the spin quantum 
Hall transition (an analogue of the quantum Hall transi- 
tion, but in symmetry class C rather than A) to the loop 
model 

For loops with crossings (p > 0) we can again con- 
struct a mapping of the first kind by taking a 'classical' 
limit in a network model with a Kramers doublet on each 
edge (replacing each quantum node with one of the three 
classical possibilities in Fig. [2]). However this correspon- 
dence is rather trivial because of the explicit suppression 
of quantum tunnelling. It would be interesting to know 
whether the analogy with localisation goes beyond this 
- in particular, whether the critical behaviour of the 
loop model can be related to the critical behaviour of a 
true localisation problem. (It is interesting to note that 
our value of v is close to estimates of v for the symplectic 
class [57].) 

Returning to loop models in their own right, there is 



a good understanding of the zoology of critical points in 
loop models without crossings, many of which fit into the 
one-parameter family of universality classes in SLE K . In 
general, crossings take us outside this family. Here we 
have discussed a new line of critical points exemplifying 
this, but we certainly do not expect that this exhausts the 
possibilities for new critical behaviour — much remains 
to be learned. 
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